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We unveil the generation of universal morphologies of fluid interfaces by radiation pressure what- 
ever is the nature of the wave, acoustic or optical. Experimental observations reveal interface de- 
formations endowed with step-like features that are shown to result from the interplay between the 
wave propagation and the shape of the interface. The results are supported by numerical simulations 
and a quantitative interpretation based on the waveguiding properties of the field is provided. 

PACS numbers: 43.25.Qp, 42.50.Wk. 



More than one century after the pioneering works of 
Poynting [1] and Rayleigh [2] on the mechanical effects 
of electromagnetic and acoustic waves, it is a common 
knowledge that both light and sound exert radiation pres- 
sure on matter. Experimentally, a striking demonstra- 
tion of its existence is the observation of deformations 
of fluid interfaces as shown in acoustics by Hertz and 
Mende in 1939 [3] and in optics by Ashkin and Dziedzic 
in 1973 [4]. Since these observations, radiation pressure 
of waves have been exploited in various contexts such as 
contactless metrology of fluids ([5] in optics and [6] in 
acoustics), liquid droplets ejection ([7] in optics and [8] 
in acoustics), and biomedical applications ([9] in optics 
and [10] in acoustics). 

Here we report on universal kind of steady state de- 
formations of initially flat fluid interfaces induced by ra- 
diation pressure whatever is the nature of the wave, be 
it acoustic or optical. We show that this novel family 
of step-like universal interface morphologies results from 
the balance between buoyancy, capillarity and radiation 
pressure for both the acoustic and electromagnetic cases. 
Performing numerical simulations we show that refrac- 
tion drives the interplay between the shape of a deformed 
interface and the wave propagation along it. Finally, we 
propose an interpretation of the observed morphologies 
based on the waveguiding properties of the deformations. 

This is made possible by using interfaces between sim- 
ple fluids for acoustic waves whereas extremely soft inter- 
faces of near critical fluids are used in the case of optical 
waves. At rest, two immiscible fluids define a planar in- 
terface, see Fig. 1(a). The interface is deformed by the 
acoustic or optical radiation pressure of a focused sound 
or light beam that impinges at normal incidence from the 
medium with the largest phase velocity (c+) to the one 
with the lowest phase velocity (c_), thereby leading to 
step-like deformations as sketched in Fig. 1(b). 

Experimentally, the acoustic beam is obtained from 
an immersed, spherical, piezoelectric transducer with a 
38.4 mm radius of curvature, 38 mm diameter operating 
at 2.25 MHz central frequency. In the optical case, we use 
a focused Gaussian beam at 514.5 nm wavelength. For 



each experiment, we set the incident power to a value 
that gives an aspect ratio H/wq of a few units, where h 
is the height of the fluid interface deformation and wq is 
the beam waist [11]. Typical examples of the observed 
morphologies are displayed in Fig. 1(c) — 1(h). In addi- 
tion, the relevant characteristics of the fluids used in the 
acoustic case are summarized in Table I whereas in optics 
we use the transparent, isotropic and nonmagnetic two- 
phase microemulsion described in [12]. At T — Tc = 4 K, 
where Tc is the critical temperature above which two dis- 
tinct phases coexist, c+ = 0.6855c and c_ = 0.6807c 
where c is the speed of light in vacuum. 

Hereafter, the observed fluid interface morphologies 
are labeled following their number n of steps, namely 
AMn and OMn in the acoustic and optical case, respec- 
tively. Quantitatively, we demonstrate that the steady- 
state morphologies obey the following equilibrium equa- 
tion between radiation and gravito- capillary pressures. 



Hr = Hgc , 
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FIG. 1. (a) Planar fluid interface at rest, (b) Step-like 
interface deformed by radiation pressure. (c,d) First (OMl) 
and second (0M2) order optical morphologies, which corre- 
spond to input beam power 470mW and 1020mW, respec- 
tively, (e-h) First (AMI), second (AM2), third (AM3) and 
higher order (AMn) acoustic morphologies, which correspond 
to input beam power 4.0W, LOW, 220mW and 20mW, re- 
spectively. 
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AMI AM2 AM3 AMn 

Fluid Water Water Salted water 25wt.% Silicone oil lOOcSt 

Fluid Kerosene Silicone oil lOOcSt Chloroform Oil FC72 

c+ (m • s-^) 1490 1490 1783 1000 

c- (m • s-^) 1315 1000 1000 512 

y9+ (kg • m-^) 998 998 1189 966 

p- (kg • m-^) 790 966 1500 1680 



TABLE 1. Fluid characteristics at room temperature for the acoustic experiments , where c± and p± refer to the phase velocity 
and the density of the fluids labeled "+" and "— " , respectively. 



In Eq. (1) the radiation pressure is = n • 

[((T+) — (T_)) • n] with n the unit vector normal to the 
interface oriented from fluid "— " to fluid and {T±) 
refers to the time averaged values of the radiation ten- 
sor at both sides of the interface over a wave cycle. 
The gravito-capillary contribution Hgc = |p+ — P-lgh — 
^£[rh\r)/^/lTh^] [13] where p± is the density of 
fluid =b, ^ is the modulus of the gravitational accelera- 
tion, 7 is the interfacial tension, and h{r) is the height of 
the axisymmetric deformed interface, with h = at rest 
and h > when the wave is turned on, and h'{r) = 
In addition, we have introduced the cylindrical coordi- 
nate system {r^(j)^z) with orthonormal basis (e^,e0,e;2), 
z being the symmetry axis of the deformation oriented 
toward its tip (see Fig. 1). 

On the one hand IIgc is calculated from the exper- 
imental interface profile h{r). On the other hand, the 
evaluation of IIr requires the computation of the field. 
This is done by solving the propagation equation for a 
given experimental profile h{r) accounting for the conti- 
nuity relations for the field at the interface and appropri- 
ate source distribution that defines the incident beam. 

In acoustics, the linear propagation equation in per- 
fect fluids expresses as A^-t — ^^^^ = where ^± 

is the complex axisymmetric acoustic pressure fleld in 
each fluid, p±{r,z,t), with time t. Assuming harmonic 
flelds with angular frequency a;, the acoustic velocity 
fleld is obtained from the Euler equation that gives 
u-t(r, = j^j^^p±- Numerically, the propagation 
problem — Helmholtz equation in a two-phase axisym- 
metric conflguration — is solved by using a boundary ele- 
ment method [14]. Once the pressure and velocity flelds 
are calculated, see panels (a) and (b) in Figs. 2, 3 and 4 
where the real part and modulus of p are shown for the 
AMI, AM2 and AM3 situations, IIr is computed using 
Brillouin expression for the acoustic radiation tensor Tac 

Tac = -^XP^l + \puH -pu0u, (2) 

where x± = {p±c±)~^ is the isentropic compressibility of 
each fluid, I is the identity tensor, and is the dyadic 
product. The validity of Eq. (1) is tested by comparing 
the radial dependence of IIr and IIgc for various mor- 
phologies as shown in panel (c) of Figs. 2, 3 and 4 [15]. 
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TABLE 11. Table of equivalence between acoustic and elec- 
tromagnetic wave propagation. 

We conclude to an overall satisfying validation of Eq. (1). 
Noticeably, the main extrema of IIr correspond to the 
ones of IIgc, which sign he step-like nature of the mor- 
phologies. We note that the observed on- axis discrepancy 
between IIr and IIgc could be ascribed to the hydrody- 
namic stress exerted at the tip of the deformation by the 
rectifled fluid flow that results from thermoviscous dissi- 
pation ("acoustic streaming") [16]. 

In the optical case, since the refractive index contrast 
is small, < 10~^, the scalar propagation equa- 

tion is known to accurately describe the electromagnetic 
fleld whatever its polarization state. In order to bene- 
flt from the mathematical equivalence between acoustic 
and electromagnetic wave propagation [17], we restrict 
our study to the particular cases of a TE-polarized fleld, 
^ = E'te where Ete = ^te^^^, or TM-polarized fleld, 
^ = i^TM where Htm = ^TMe^/,, where E and H are 
the electric and magnetic flelds. This allows to use the 
numerical code developed for the acoustic case following 
the polarization dependent correspondance summarized 
in Table II. Once the electric and magnetic flelds are cal- 
culated, see panels (e) and (f) in Figs. 2, 3 and 4 that 
show the real part and modulus of £^te for the OMl and 
0M2 situations, IIr is computed using Maxwell expres- 
sion for the electromagnetic radiation tensor Tem [18] 

Tem = -^^E^I - ^/^H^^ + eE E + /iH (8) H , (3) 

where e± = eo \fcfc± and pL± = /io are the dielectric and 
magnetic susceptibilities of each fluid, respectively, /io 
and eo being the ones of vacuum. From the panels (g) and 
(h) of Figs. 2 and 3, the comparison of IIr and IIgc [15] 
allows to draw the same conclusions as in the acoustic 
case, including the possible role of a light induced on- 
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axis flow ascribed to light scattering by refractive index 
critical fluctuations ("optical streaming") [19, 20]. As 
shown in Fig. 2(e), where IIr is plotted for both TE and 
TM cases, the radiation pressure is almost insensitive to 
the polarization state of the wave, as expected. 

In previous works we demonstrated that the waveguid- 
ing of the fleld along the deformed interface constitutes 
the feedback mechanism for radiation pressure effects. 
In particular, this can lead to the stabilization of trans- 
lationally invariant (here along the z axis) liquid columns 
by acoustic [21] or electromagnetic [22] guided waves. 
Next, we show that this mechanism applies to noninvari- 
ant fluid deformations as well. We note that 0M2 and 
AM2 deformations have been already observed, however 
not explained, in [23] and [24], respectively. 

For this purpose, we compare (i) the numerically com- 
puted fleld at altitudes that correspond to locally cylin- 
drical shape of a given deformation with radius R with 
(ii) the fleld of the dominant guided mode existing for 
a cylindrical waveguide with same radius. As shown 
in Figs. 5, 6 and 7, we observe a striking coincidence 
that leads us to calculate the modal content of the fleld 
propagating along the deformation. The 2;-dependence 
of the modulus of the normalized amplitude \an\ of the 
nth guided mode [25] is shown in Figs. 5(a) and 5(d), 
6(a) and 6(e), and 7(a) for the AMI, AM2, AM3, OMl 
and 0M2 situations. Note that the fundamental mode 
n = 1 of a cylindrical waveguide exists whatever its ra- 
dius whereas higher order modes n > 1 exist above cut-off 
radii r = Rn-i, respectively. The observed plateaus of 
\an\ versus z coincide with locally z-invariant morpholo- 
gies. This indicates that every cylindrical portion of the 
fluid deformation actually behaves as a waveguide. 

Interestingly, the step-like morphologies can be clas- 
sifled using a dimensionless quantity that only depends 
on the wavelength, the beam waist and the phase ve- 
locities of each fluids. Following standard waveguiding 
theory [26], we introduce the characteristic normalized 
frequency Vc = A:+i(;o\/(c+/c_)^ — 1, where A:+ is the 
wavevector in fluid . The number of existing guided 
modes is known to increase with Vc, so does the order 
of the observed morphologies. Indeed Vc = 5.0, 6.0, 8.1, 
and 9.0 for the acoustic morphologies AMI, AM2, AM3 
and AMn shown in Fig. l(e)-(h), respectively, whereas 
Vc = 3.6 and 6.1 for the optical morphologies OMl and 
0M2 shown in Fig. 1(c) and 1(d), respectively. However, 
we notice that a given morphology guides a larger num- 
ber of modes in optics than in acoustics. This can be 
grasped by comparing, on the one hand, the cases AMI 
and OMl, and, on the other hand, the cases AM2 and 
0M2. This could be ascribed to different values of kwo, 
which is a constant of the order of 1 in acoustics whereas 
it is typically one order of magnitude larger in optics. 

We conclude to the observation of universal step-like 
morphologies of fluid interfaces deformed by the radia- 
tion pressure of acoustic and electromagnetic waves. This 



phenomenon basically relies on the waveguiding proper- 
ties of sound or light in axisymmetric two-phase liquid 
deformations with distinct phase velocities. 
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FIG. 2. (color online) Analysis of the AMI and OMl cases, see Fig. 1(e) and Fig. 1(c). In each case, the spatial distribution 
of the real part and the modulus of the field in an equatorial plane are shown with the deformed interface profile superimposed 
on it as a white solid curve [panels (a,b) and (e,f)]. (c,g) Balance between the radial dependence of the radiation pressure (IIr) 
and the gravito-capillary contribution (IIgc), see [15] for details. (d,h) Deformed interface profiles. The altitude at the tip of 
the deformation is set here to zero in all cases. 
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FIG. 4. (color online) As in Fig. 2 for the AM3 case, see 
Fig. 1(g), and [15] for details. 



6 



Acoustic Morphology 1 Optical Morphology 1 




FIG. 5. (color online) Analysis of the AMI and OMl cases, see Fig. 1(e) and Fig. 1(c). (a,d) ^-dependence of the modulus 
of the normalized amplitude \an\ of the nth guided mode. (b,e) Deformed interface profiles. The altitude at the tip of the 
deformation being set to zero in all cases. Each colored areas refers to radial range of existence of a given set of guided modes. 
(c,f) Correspondence between the numerically computed field at altitudes that correspond to locally cylindrical portion of a 
given deformation (solid curve) with the field of the dominant guided mode existing for a cylindrical waveguide with same 
radius (dashed curve), with ip{r, z,t) = ^ ^ ' ^LT^ (r^z^'t)!^ ' ^ ' where ip = p in acoustics and ip = ^te in optics. 





